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A dynamical appearance of scaling solutions in the relativistic hydrodynamics applied to describe 
ultra-relativistic heavy-ion collisions is studied. We consider the boost-invariant cylindrically sym- 
metric systems and the effects of the phase transition are taken into account by using a temperature 
dependent sound velocity inferred from the lattice simulations of QCD. We find that the transverse 
flow acquires the scaling form r/t within the short evolution times, 10-15 fm, only if the initial trans- 
verse flow originating from the pre-equilibrium collective behavior is present at the initial stage of 
(-^ ■ the hydrodynamic evolution. The amount of such pre-equilibrium flow is correlated with the initial 

pressure gradient; larger initial gradients require smaller initial flow. The results of the numerical 
calculations support the phenomenological parameterizations used in the Blast- Wave, Buda-Lund, 
and Cracow models of the freeze-out process. 
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^ . T. INTRODUCTION 

y—i 

Hadronic data collected in the heavy-ion experiments at RHIC support the idea that the system formed in Au+Au 
collisions is highly thermalized and undergoes strong transverse and longitudinal expansion Moreover, successful 
parameterizations of the freeze-out process indicate that such expansion may be characterized by the Hubble law 
(2, 0,0] I n the simplest form, used in cosmology, this law states the proportionality of the relative velocity of galaxies 
to their relative separation, 

§ I v = Hr. (1) 

The constant of proportionality or Hubble's constant H is a function of time. In the Friedmann universe [f| as well 
I as in analytic solutions of non-relativistic fireball hydrodynamics @ , its value is 

^ : h=i 

where the function R(t) is in general a complicated function of time. In nuclear hydrodynamics, the Hubble law 
characterizes the fluid velocity distributions of the expanding matter, while in cosmology the Hubble law characterizes 
the expansion of the space. In the non-relativistic hydrodynamics, the scale parameter R(t) depends on the initial 
conditions as well as on the equation of state, while in cosmology the time evolution of the scale parameter R depends 
not only on the initial conditions (flat, open or closed universe) and the properties of matter (matter dominated 
and radiation dominated epoch) but also on the possible existence of dark energy and cosmological constants, and 
the possible presence of an exponentially accelerating, inflatory period. At the end of the accelerating period, when 
R = const, and R « Rt, the Hubble constant is simply the inverse of the lifetime, H = l/t. 

At the very end of fireball explosions, the pressure decreases to vanishing values, hence the acceleration caused by 
pressure gradients becomes negligibly small. In this case, the Hubble law connects the hydrodynamic flow four- velocity 
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with the space-time position of the fluid element s* 1 , in a simple way, 




(3) 



The quantity r in Eq. is the proper time characterizing the freeze-out hypersurface, 



t = \J t 2 — r 2 — z 2 , r = V 1 x 2 + y 2 . 



(4) 



Clearly, the parameterization J3J - J3} makes sense in the space-time region defined by the condition r 2 + z 2 < t 2 . 
The three-velocity field of the form 



following directly from Eq. , appears very often in the analysis of the hydrodynamic equations applied to describe 
hadronic or nuclear collisions. In this context it is called the asymptotic scaling solution 0, IE Hi E3> El 0] In 
particular, for the boost-invariant systems the l ong itudinal velocity must be of the form v z — z/t, which is a direct 
consequence of the Lorentz symmetry [Illll3lll4lll5lll6lll7l |. It is often believed, that during the pre-equilibrium period 
of high energy nuclear collisions, no significant transverse flow is gene rated (for recent reviews of the hydrodynamic 
description of relativistic heavy- ion collisions see, e.g., Refs. [Illll El)- In such a case, the transverse fluid 
velocity builds up during the hydrodynamic evolution of the system ll '2f and the scaling form v r = r / 1 may be reached 
only for sufficiently large times (with details depending on the equation of state and initial conditions). Similar features 
characterize also a spherically symmetric expansion of the system being initially at rest. In this case, the radial flow 
is formed by the outward action of the pressure gradient, and the numerical calculations show that Eq. J3J is the 
asymptotic solution 0, . 

Recently, there is a revived theoretical interest in finding exact solutions of relativistic hydrodynamics. Biro found 
an interesting solution, relevant for the case of a vanishing speed of sound, which interpolates between an early Bjorken 
type of the flow profile and the Hubble profile in the late period of the expansion |22l |23| . New exact solutions of 
relativistic hydrodynamics were found, using more general equations of state, in the 1+1 dimensional case |24j and 
in the 1+3 dimensional case; for axially symmetric H^,|2^|, as well as for ellipsoidally symmetric expansions |27| . 
Although these solutions contain arbitrary scaling functions describing the rapidity profile, the flow profile in all of 
these works coincides with the Hubble law. 

In view of the success of the fits 0, 0, |2^, HiJ which all indicate very strong transverse flows, one of the central 
issues is whether the times actually available in relativistic heavy-ion collisions are sufficient to allow for a dynamical 
development of such strong transverse flows, in particular, the scaling solutions corresponding to Eq. Q. The 
situation is especially intriguing for RHIC, where several measurements indicate an unexpected, rather short (about 
10 - 15 fm) duration time of the collision process. For example, one finds r ~ 10 fm using the RHIC data in the 
relation ^(mj) = T^Tk/mr |2lJ, which connects the longitudinal pion correlation radius Rl, the kinetic freeze-out 
temperature the evolution time r, and the transverse mass of the pion pair mx- 

These measurements, when interpreted with care, yield only the inverse of the (longitudinal) Hubble constant, which 
can be identified with the lifetime only if the scaling solution is assumed to be developed in the longitudinal direction. 
This situation is analogous to the estimate of the lifetime of the Universe: the inverse of the presently measured 
value of the Hubble constant yields an order-of-magnitude estimate of the lifetime, which has to be corrected for the 
effects of inflation and possible other acceleration periods. Similarly, in nuclear hydrodynamics, there is an initial 
longitudinal acceleration period, hence the estimated lifetimes should be interpreted only as (lower) limits on the total 
lifetime of the system 

Many hydrodynamic codes used to describe the heavy-ion data show that the scaling solutions do not appear before 
the freeze-out of the system. However, such approaches assume commonly that the initial transverse flow is zero. An 
exception from this rule is the work by Kolb and Rapp |3l| . where the pre-equilibrium transverse flow is considered 
and its presence improves the agreement of the model calculations with the data. Another important exception is 
a class of the non-relativistic, selfsimilar solutions of the fireball hydrodynamics, which is by now solved completely 
in the ellipsoidally symmetric case for abritrary initial sizes and expansion velocities of the principal axes of the 
expanding ellipsoids |32j|. arbitrary initial temperature profiles |33f . as well as for arbitrary (temperature dependent) 
speed of sound [34| . 

In this paper we follow such ideas and assume that the elementary parton-parton collisions, leading to the thermal- 
ization of the system, lead also to collective behavior and development of the transverse flow already at the initial stage 
of the equilibrium hydrodynamic evolution at the time t = to ~ 1 fm. For simplicity, we consider the boost-invariant 
and cylindrically symmetric systems with the initial transverse flow defined by the formula 




(5) 



v r = v r (t =t a ,r) 



Hr 



(6) 



Vl + H 2 r 2 ' 
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The parameter H in Eq. 10 may be interpreted as the Hubble constant which determines the magnitude of the 
initial transverse flow; in the range r < 1/H the flow is well approximated by the linear function v® ~ Hr, whereas 
for r > 1/H the flow approaches the speed of light, a boundary condition frequently assumed in the hydrodynamic 
equations for large values of r 12]. Also, if the above equation is specified on a hypersurface corresponding to a 
constant proper-time, it coincides with the Hubble law of Eq. @. 

By varying the value of H we control the amount of the initial transverse flow which may possibly develop into the 
scaling form JSJ). We note that in the scaling solution (JjjJ the role of the Hubble constant is played by the inverse of 
the time coordinate t. This means that setting H = 1/to = lfm -1 we assume that the initial flow is already very 
close to the scaling form at to- In this particular case the question arises if the flow remains close to the scaling form 
in the subsequent time evolution of the system. 

A few comments are in order now. Following the conventional hydrodynamic approach |l2j . we specify the initial 
conditions at a constant time, to = 1 fm, and search for the solutions which are regular functions of t and r in the 
region: t > to, r > 0. In this case, at a given value of time t, the scaling solution cannot hold for arbitrary large 
values of r, since this would yield the flow velocities larger than the speed of light, and also the initial condition 
deviates from Eq. (1) for large values of the transverse radius. Hence, in our analysis we search for the solutions of 
the hydrodynamic equations which yield the flow profiles possibly close to the scaling solution in the region r < t. In 
a separate paper we intend to explore a different type of the initial conditions, which are specified at a fixed value of 
the proper time and lead, in certain special cases, to the exact solutions given in Refs. |2(|[27|]. We also note that in 
view of the recent BRAHMS data HI HI describing rapidity dependence of hadron production, the boost-invariant 
approach assumed in the present calculation may be appropriate only if limited to the rapidity range: — 1 < y < 1. 



II. CHARACTERISTIC FORM OF THE HYDRODYNAMIC EQUATIONS 

In this Section we introduce the basic notation and rewrite the hydrodynamic equations in the form convenient for 
the numerical calculations. We follow closely the method introduced by Baym, Friman, Blaizot, Soyeur, and Czyz 
|12|. We restrict our considerations to the systems with zero net baryon density, which is a good approximation for 
description of the central rapidity region at the RHIC energies (thermal analysis of the ratios of hadron multiplicities 
indicates that the baryon chemical potential at RHIC energies is about 30 MeV, i.e., it is much smaller than the 
corresponding temperature of about 170 MeV |4l 1371 138| ) . In this case the hydrodynamic equations have the form 

u»d»{Tu v )=d v T, (7) 



9 M (au») = 0, (8) 

where T is the temperature, a is the entropy density, and = 7(l,v) is the hydrodynamic four- velocity (with 
7 = 1 /yl — v 2 ). Eq. (7J is the acceleration equation which is the analog of the Euler equation of the classical 
hydrodynamics, whereas Eq. JHJ represents entropy conservation (the adiabaticity of the flow). Since T is the only 
independent thermodynamic variable, all other thermodynamic quantities can be obtained from the equation of state 
P (T), connecting pressure P with the temperature T. With the help of the thermodynamic relations 

de = Tdcr, dP = adT, w = e + P = T<r, (9) 

other thermodynamic quantities, such as the energy density e or the enthalpy density w, can be obtained. In addition, 
the equation of state allows us to calculate the sound velocity 

c '=aF = T^- (10) 

Equation (jHJ and the spatial components of Eq. (JJJ) may be rewritten for the boost invariant systems with cylindrical 
symmetry as 

<91nT d\nT da da 
Vr—^ + —^ \--Z7+ v rs- = 0, (11) 



dt dr dt d 



r 



91ner dlna da da 1 v r 

dt 7 dr r dt dr t r 

where a is the transverse rapidity of the fluid defined by the condition v r = tanha. The longitudinal component has 
the well known boost- invariant form v z = z/t |ll| . 
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FIG. 1: a) Temperature dependence of the sound velocity as obtained from different theoretical models described in the text, 
b) Sound velocity used in this calculation: the solid line describes the lattice result, i.e., the function "lattice I" from part a) 
extrapolated to low temperatures, whereas the dashed line describes the analytic model defined by Eq. 12ft . 



By introducing the potential $ (T) defined by the differentials 



rflnT 

d<l> = = c s a In er, 



(13) 



and by the use of the two functions a± defined by the formula 

a± — exp ($ ± a) , 
Eqs. Ijlll) and (|12|l may be cast into the characteristic form ^1 

— a ± (t r) + Vr±Cs d a± U r ) + Cs 
dt ' 1 ± v r c s dr ' 1 ± v r c.s 



- + 7)«± (*,r)=0. 



If the functions a± are known, the potential $ may be calculated from the formula 



$ = iln (a + a_) 



(14) 



(15) 



(16) 



and the velocity is obtained from the equation 



a + + a— 



(17) 



The knowledge of the function c s (T) allows us, by the integration of Eq. (|13f) . to determine $ as a function of the 
temperature; this function will be called later <!>t- However, to get a closed system of equations for a + and a_, we 
have to invert this relation and obtain T as a function of <&; this function will be called later T$. In this way, the 
sound velocity may be expressed in terms of the functions a + and a_, 



T, v . ( iln (a . a ) 



(18) 



and Eqs. (|15|) may be solved numerically. 



III. TEMPERATURE DEPENDENT SOUND VELOCITY 



In our numerical calculations we take into account the temperature dependence of the sound velocity. In this way, 
we generalize the approach of Ref. where Eqs. Ijl5|l - l|17(l were solved numerically only in the case c 2 s — 1/3. 

The results of different model calculations of the sound velocity, which may serve as the input for the hydrodynamic 
calculations are presented in Fig. 2] a). The solid line (denoted as "lattice I") is the result of Mahonty and Alam |39| . 
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who compiled the lattice results obtained by Karsch |4(| in order to get c s {T) from the temperature dependence of 
the energy density. The long-dashed line (denoted as "lattice II") shows the result of the lattice QCD calculations by 
Szabo and Toth g^. 

A sudden but smooth change of the sound velocity in the small temperature range around T = T c , as observed 
in the lattice calculations, see Fig. n a )> indicates a rapid cross-over from a hadron gas to a quark-gluon plasma 
phase. Above the critical temperature, (T 3> T c ) the sound velocity approaches the limit valid for massless particles, 
c 2 s = 1/3, whereas below the phase transition (T <C T c ) the sound velocity is much smaller, being close to the value 
obtained for the case of the noninteracting gas of hadron resonances. 

The hadron-gas result shown in Fig. ED a) was obtained by us in the calculation which uses the recent fits to the 
meson and baryon mass spectra |42l I43.l44l| denoted below by pm(iti) and pB{m). One possible parameterization of 
the spectra 0], which reveals directly their exponential growth, as suggested long ago by Hagedorn (4^, is as follows 



PM(m) = a M exp(77i/T M ), a M = 4.41 GeV \ T M = 311 MeV, 

p B (m) = a B cxp(m/T B ), a B = 0.11 GeV -1 , T B = 186MeV. (19) 

With the help of the mass distributions (|19fl . we calculate the entropy density of the hadron gas as a sum over 
contributions from all known hadronic states from the formula 



a{T) = 2^J pM(m)m 3 K 3 (-J dm+—j p B (m)m 3 K 3 (-J dm. (20) 

The sound velocity of the hadron gas follows then directly from the use of Eq. (|2L)fl in I|1U|) . Eq. I|2U|) is valid in the 
case of zero baryon chemical potential and the factor 2 in the second term accounts for antibaryons. The effects of 
the quantum statistics (Bose-Einstein or Fermi-Dirac) are neglected in Eqs. Q20[l . because they are known to be small; 
of the order of 20% or less for reasonable range of the temperatures @ . To match our hadron-gas calculation with 
the lattice data we assumed that the critical temperature is 170 MeV. The upper limits of the integrations in H20|) are 
-^mesons = 2.3 GeV for mesons and Mg^* ons — 1.8 GeV for baryons. These limits are determined by the range where 
the fit l|19|l works well |4J|; for higher masses, due to the lack of data, the spectra saturate. 

For comparison, in Fig. ^ a ) we show the speed of sound obtained in the similar calculation where the massive 
pions are included only. The speed of sound in the pion gas reaches very fast the limiting value of 1/ V3, which may 
be confronted with the non-monotonic behavior of c s in the gas of resonances. Similar behavior of the speed of sound 
in the gas of resonances was found also in the case with non-zero baryon density 0] . It is interesting to note that the 
lattice data agree with the hadron-gas calculation close to the phase transition region if we assume T c ~ 170 MeV. 
Moreover, the speed of sound in the hadron resonance gas below the critical temperature is found to be significantly 
smaller than which is the limit of massless ideal relativistic gas used in the bag model type of the equations 

of state. As the hydrodynamic equations describe the properties of matter through the equation of state, or more 
precisely through the speed of sound @, such decrease of the speed of sound in the hadron gas stage, compared to 
the massless pion gas limit, changes drastically the corresponding time evolution of the hydrodynamic solutions. 

Although the calculations of the speed of sound shown in Fig. are based on the observed hadron mass spectra 
below T c and on the lattice QCD above T c , they are still somewhat ambiguous and not completely satisfactory. 
For example, lattice QCD calculations below the critical temperature still yield too heavy pions, and we know that 
the value of the speed of sound is rather sensitive to the pion mass. On the other hand, the calculation based on 
the Hagedorn mass spectrum in the hadronic phase is more realistic than the lattice results in the low temperature 
domain, however, it neglects the role of interactions between the hadrons. 

In this situation, we decided to use as the main input the equation of state (sound velocity) delivered by Ref. [3^, 
which is shown as "lattice I" in Fig. [I] a). This calculation is based on first principles and extends down to about 
0.6 T c , i.e., to the low-temperature region of about 100 MeV relevant for the discussions of the kinetic freeze-out. In 
the actual calculation we extrapolate this result to even lower temperatures as shown in Fig. \}\b). We also use an 
analytic form of the function c s (T) which exhibits the main features observed in Fig. Q]a) and, at the same time, 
leads to the analytic expressions for the functions &t(T) and T$($). This function is defined by the formula 



i-y 1 



(21) 



2 V 1 + (T/T) 2n , 

Using Eq. (J3TJ one finds c s (T) = 1/V3 for T > f and c s (T) = l/(2\/3) for T«f. A straightforward integration of 
Eq. (|13|) gives in this case 

M T) = ^ln V' f r (22) 
2n 1 + 2(T/T) 2 " 
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FIG. 2: a) The potential $ obtained by the integration of Eq. 113H . b) Temperature as a function of the potential In both 
cases, a) and b), the solid lines show the results obtained with the lattice equation of state, whereas the dashed lines describe 
the analytic model defined by Eq. (1211 . 



and 



e ^ 



:7i<l> 4n<f> 

e v-3 -|- e \/3 



(23) 



With the parameters T = 1.08 T c and n = 6 the function Q21JI well approximates the lattice results I and II in 
the region slightly above T c , and interpolates between the lattice results I and II at higher temperatures. At low 
temperatures, T < 0.6 T c , the function l|21l) behaves like a constant, see Fig. \I\b). The functions $t(T) and T$($), 
defined by Eqs. Q22J1 and l)23|l with T — 1.08 T c and n = 6 are represented in Fig. [3 by the dashed lines. The solid 
lines in Fig. [3 show the same functions obtained for the lattice case. (In the following, we shall refer to the lattice 
calculations having in mind the case "lattice I", including a linear extrapolation at very low temperatures). 



IV. INITIAL CONDITIONS 



From the symmetry reasons, the velocity field should vanish at r = 0. This condition is achieved if the functions 
a+ and a~ are initially determined by a single function a(r) according to the prescription |T^ 

a+(t = t ,r) = a(r), a_(t = t ,r) = o(-r). (24) 

In the following we shall assume that the hydrodynamic evolution starts at a typical time t = to = 1 fm. We shall 
also assume that the initial temperature profile is connected with the nucleon- nucleus thickness function Ta(t), 

T(t = t ,r) = const T\ /3 (r). (25) 

where 

T A (r) = 2 J d 2 s J dzp(V(s-r) 2 +z 2 ) . (26) 

Here the function p(r) is the nuclear density profile given by the Woods-Saxon function with a conventional choice of 
the parameters (p = 0.17 fm~ 3 , r = (1.12A 1 / 3 - 0.86A" 1 / 3 ) fm, a = 0.54 fm, A = 197). 

The idea to use Eq. I|25|l follows from the assumption that the initially produced entropy density o ~(r) is proportional 
to the number of the nucleons participating in the collision at a distance r from the collision center 19], cr(r) ~ Ta(t). 
Since for massless particles the entropy density is proportional to the third power of the temperature, we arrive at 
Eq. (25). 

We note that other choices for the initial temperature profile are also conceivable. If we assume that the initially 
produced energy density is proportional to the nuclear thickness function, instead of Eq. I|25|l we obtain 



T(t = t , r) = const T^ /4 (r). 



(27) 
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FIG. 3: Hydrodynamic expansion of matter being initially at rest, i.e., in the case H = 0. The initial central temperature 
To = 2T C . The part a) shows the transverse velocity as a function of the distance from the center for 6 different values of time, 
ti = 1, 4, 7, 10, 13 and 16 fm. The dashed thin lines describe the ideal Hubble-like profiles of the form r/ti (r < ti). The part b) 
shows the transverse four-velocity as a function of r for the same values of time. The part c) shows the temperature profiles 
in r whereas the part d) shows the isotherms in the t — r space. In this case, the labels at the curves denote the values of the 
temperature in units of the critical temperature. 



If we assume that the energy deposition into the thermalization is driven by the collisions of wounded nucleons, and 
every collision contributes with certain probability distribution to a local increase of the temperature, then after many 
collisions the central limit theorem may describe the form of the local temperature distribution and this is a Gaussian 
form. However, if there are big fluctuations in the deposited energy, the generalized central limit theorems apply and 
in this case the local temperature distribution may have the generalized, Levy stable form. As a special case, the 
Lorentzian temperature profile can also be obtained. 

The two initial conditions JSJl and 1|25|) may be included in the initial form of the function a(r) if we define 



a(t = t ,r) = a T (r) 



(28) 



where 



a T (r) = exp $ T (const T|(|r|) 



(29) 



V. RESULTS 



It is important to observe, that the Hubble law v = Hr is satisfied in the r <C t region after t = 7 fm in all the 
cases that we explored numerically in the present calculation, regardless of the initial conditions. However, the value 
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FIG. 4: Hydrodynamic expansion of matter with initial pre-equilibrium flow characterized by the velocity profile © with 
H = 0.25. The initial central temperature To = 2T C . Notation as in Fig. |3] 

of the Hubble constant typically deviates form the inverse of time, which signals that the asymptotic form of the 
Hubble flow, eq. is reached only after a longer time period. The onset of the asymptotic scaling is determined by 
the initial conditions as we shall detail below. 

In Figs. |21 and 01 we show our results obtained for two different initial conditions characterized by the two different 
values of the parameter H, H = and H — 0.25 fm -1 , respectively. In both cases the lattice equation of state is used 
and the initial temperature in the center of the system is assumed to be twice the temperature of the phase transition, 
To = Tit = to,r = 0) = 2T C . This means, that for the commonly accepted value of T c , which is about 170 MeV, the 
initial temperature in the center reaches about 340 MeV. 

In the case H = 0, the transverse flow is initially set to zero but it builds up during the evolution of the system, 
as shown in Fig. |3a). The corresponding values of the transverse four- velocity u r are plotted in the part b). To 
check if the flow approaches the asymptotic scaling solution, we compare the velocity profiles calculated numerically 
at different times U = 1,4,7,10,13,16 fm (solid lines) with the ideal Hubble-like velocities of the form r/U (thin 
dashed lines) in the regions r < U. As expected, in the case H = the calculated profiles do not agree with the ideal 
curves in the considered evolution times. 

In the part c) of Fig. |3] we observe that the central part of the system cools down very fast from T = 2T C down 
to the critical temperature T = T c and the subsequent cooling is very much slowed down. This behavior is caused 
by different values of the sound velocity in the regions above and below T c ; larger values of c s in the plasma phase 
imply its faster cooling. We note that for the first order phase transition the speed of sound drops to zero at T = T c 
and the system keeps on expanding at a constant temperature of T — T c . In the present case we deal with a sudden 
cross-over transition, hence the expansion of the volume is coupled to a small decrease of the temperature. From 
Fig. |21 d) showing the isotherms in the t — r space, we may conclude that expansion of the system without noticeable 
cooling below T c takes more than 20 fm. 

The evolution of matter shown in Fig. |2]may be confronted with the situation where the non-zero pre-equlibrium 
flow is included in the initial condition. In Fig. 01 we show our results obtained in the case H = 0.25 fm" 1 , with 
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the same initial central temperature To = 2T C . Since the transverse flow is present already at the beginning of the 
evolution, the expansion of the system is much faster than that discussed in the previous case of H = 0. In Fig. 21a) 
we show the velocity profiles in r, again for 6 different values of time. By comparing to the ideal curves of the form 
r/U (thin dashed lines) we observe that the flow approaches the asymptotic scaling solution after about 7 fm. 

In Fig. |2| c) we can see similar behavior to that presented in Fig. c), namely, the initial fast cooling of the hot 
center, which is slowed down when the system approaches T c . We observe, however, that in the case H = 0.25 fm - 
the slowing down of the cooling process is not as much effective as that observed in the case H = 0. Due to the larger 
transverse flow, the energy from the interior is transported outside, the temperature in the center continues to drop 
down, and a very interesting situation happens: the parts of the system away from the center become hotter than 
the parts in the center. This phenomenon is well seen in Fig. 0] d) where the isotherms of the system are shown. 
We note that isotherms of similar shape are used in the Cracow model where they are defined by the condition of 
constant proper time r • On the other hand, the Blast- Wave model assumes that freeze-out happens at a constant 
value of the ordinary time t with a fixed temperature T; whereas the Buda-Lund model assumes that the temperature 
profile may eventually decrease to at large transverse extensions, hence capturing the feature that the temperature 
vanishes for very large transverse coordinates in all of the presented calculations. 

In Fig. [5] we show a collection of the velocity profiles obtained for four different values of the parameter H and for 
three different values of the central temperature To. The four rows of smaller figures describe the results obtained with 
H = 0.0, 0.10, 0.25, 1.0 fm , while the three columns correspond to the central temperature To = 1.5T C , To = 2T C , 
and T — 3T C . 

The results presented in Fig. [SJare obtained with the analytic model for the temperature dependence of the sound 
velocity, hence, by comparison of Fig. [SI with the two previous figures, the effects of the change of the function c s (T) 
on the time evolution of the system may be observed. For example, comparing the part c2) of Fig. [SJwith the part 
a) of Fig. 0] (both results obtained with H — 0.25 fm -1 and To = 2T C ) we can see that the flow obtained with the 
analytic model is closer to the asymptotic scaling form than the flow obtained for the lattice equation of state. We 
note that the main difference between the two cases is that the sound velocity in the analytic model is smaller in 
the hadronic phase. The asymptotic scaling solutions are more easily generated at the reduced values of the sound 
velocity as shown in Refs. [lOl |22J, |23j . This effect explains a better agreement of the generated transverse flow with 
the asymptotic scaling solution obtained in the analytic model. 

Let us now discuss the impact of the initial temperature on the formation of the transverse flow. Since the spatial 
extension of the system is the same in all considered cases (roughly a diameter of the gold nucleus) , a higher initial 
temperature implies a larger pressure gradient leading directly to the formation of the stronger flow. This effect is 
seen in Fig. [5] if the results obtained with the same value of H but different values of To are compared. The presence 
of the pre-equilibrium transverse flow also helps to develop the strong transverse flow, and this effect is added to the 
effects of the pressure gradient. This is clearly seen in the case H = 0.10 fm -1 , depicted in the parts bl) - b3) of 
Fig. [SJ At t = 16 fm the flow in the range r < 10 fm is close but below the asymptotic scaling solution for To = 2T C , 
and close but above the asymptotic scaling solution for T = 3T C . For T = 1.5T C the flow is noticeably below the 
asymptotic scaling form. We conclude that the asymptotic solution may be reached either from above or from below 
(in the region r << i), depending on the value of the initial temperature. If the pre-equilibrium flow is strong enough 
the asymptotic solution is approached from above. This behavior is indicated by our results obtained in the cases 
H = 0.25 far 1 and H = 1.0 far 1 , Fig. EJcl) - d3). Interestingly, in the case H — l/to — 1.0 ha , where the initial 
flow agrees with the scaling form already at the beginning of the time evolution, the pressure gradient accelerates 
the matter and the convergence to the scaling form is delayed. On the other hand, the existence of a non-zero pre- 
equilibrium flow seems to be a necessary condition for the formation of the accelerationless Hubblc-like flows in the 
evolution times of 10-15 fm. This is indicated by our results obtained with H = for different values of To, Fig. 
al) - a3). The amount of the pre-equilibrium flow required to achieve the fast convergence to the asymptotic solution 
depends on the initial temperature (pressure gradient); smaller values of H may be compensated by larger values of 
T). For example, at t = 16 fm the flow profiles are very similar in the cases: H = 1 far 1 and T = 1.5Tq (part dl), 
H = 0.25 far 1 and T = 2T C (part c2), and also H = 0.1 far 1 and T = 3T C (part b3). It is quite remarkable, that 
the presence of this pre-equilibrium flow is required only to set the value of the Hubble constant to H = 1/t i.e. to 
reach the asymptotic scaling solution within a short time period, however, we find that the linear flow profile, v = Hr 
develops with H ^ 1/t in all the considered cases by about 7 fm/c, regardless of the details of the initial conditions. 

VI. CONCLUSIONS 

Our results show a dynamical development of the scaling solutions in the relativistic hydrodynamics applied to 
relativistic heavy-ion collisions. For the realistic initial conditions connecting the amount of the initially produced 
entropy with the number of participating nucleons, we find that a Hubble type linear transverse flow, characterized 



10 



a1) a2) a3) 






16 




3 








H =0.10 




T„ = 1.5T C 





;= ^^i 

16 


1 

0.8 






0.6 






> 




3 


0.4 




H =0.10 


0.2 




T„ = 2T C 












Cl> c2> c3> 


























H =1.00 






T„ = 1.5T C 


5 


10 


15 20 




r[fm| 





1 

0.8 


MM 






1 

0.8 


0.6 






0.6 


> 








> 


0.4 








0.4 


0.2 






H =1.00 


0.2 








To = 2T C 

















5 


10 


15 20 




10 
r[fm| 



FIG. 5: Velocity profiles in r (solid lines) calculated for different values of the parameter H and different values of the initial 
central temperature To. The labels 1,4,7,10,13 and 16 denote the evolution time in fm. The values of H are given in 1/fm. The 
thin dashed lines show the ideal profiles of the form r/t obtained for the same values of time. 



by v — Hr develops by about 7 fm/c regardless of the details of the initial conditions when varied numerically in a 
reasonable range. However, it is more difficult for the transverse flow to achieve the asymptotic, accelerationless scaling 
form: to reach the values of H = 1/t within the evolution time of about 10 -15 fm. For this, the necessary conditions 
are thus stronger: pre-equilibrium transverse flow has to be present already at the beginning of the hydrodynamic 
evolution at the initial time of 1 fm. The amount of the pre-equilibrium flow, required for the fast approach to the 
asymptotic solutions, depends on the initial pressure gradients; for larger gradients smaller initial flows are necessary. 

The results of these calculations give support for using the phenomenological parameterizations of the freeze-out 
process such as the Blast- Wave, Buda-Lund, and Cracow models. Certainly, more work should be done to combine 
the output of the hydrodynamic calculations with the description of freeze-out used in these models, however, we 
showed that such features as the linear flow profiles or the isotherms defined by the constant value of the proper time 
appear as the solutions of the relativistic hydrodynamics with suitably chosen initial conditions. 
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